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Abstract 

New multi-color UBVR light curves of the eclipsing binary KRCyg were 
obtained in 2005. Photometric solutions were derived using the Wilson- 
Devinney method. The result shows that KRCyg is a near-contact binary 
system with a large effective temperature difference between the components, 
approximately 5230 K. All the times of minimum light were collected and 
combined with our observations obtained in 2010 and 2011. Analysing all 
the times of mid-eclipse, the period of quasi-sinusoidal variation is about ~76 
years. The periodic oscillation could be explained by the light-time effect due 
to a presumed third component. 
Keywords: 

(stars:) binaries: eclipsing — (stars:) binaries: general — stars: early-type 
- stars: individual (KRCyg) 

1. Introduction 

Near-contact binaries are very important source in understanding the 
formation and evolution of the binary systems. The members of these sys- 
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terns have orbital periods less than 1 day. Most of the systems contain a 
component at or near its Roche lobe. Theoretical studies suggest that the 
near-contact binaries could be in th e intermediate-stage between a detached 



or semi-detached configuration (e.g. 



Hilditch et al 



1988 



Shaw 



199(1 ) . There- 



fore, these systems are important observational targets which could play a 
key role in understanding the evolutionary stages of interacting binaries. 
KRCyg (HP 333645, V = 9 m .19) was recognized as a variable star by 



Schnellerl (11931a 



cj). Its properties are relatively poorly k nown compared 



to those of other short-period binaries. The observations of 



Vetesnik! fll965f ) 



revealed the fir s t phot oelectric light curves and a period of d . 8451538 for the 



system. 



Shawl (Il990h included KRCyg in their list of near-contact eclips- 



ing binaries. Its more detailed history until 2004 has been described by 



Sipahi fc Giilmenl (120041 ) . Unfortunately, no spectroscopic study exists. 



We followed the observations of the system and obtained new light curves 
in 2005. We report the new photometric observations, and the solutions of 
the light curves. In this paper, we investigate possible implications of the 
orbital period change of the system based on our new minima times, and 
on historical data collected from the literature. Both the period analysis 
and the light-curve analysis show that KRCyg is probably a triple system. 
The system is i mportant not only because it joins the group of near-contact 
binary systems ( iShawl . Il990l ). but also because we show that it is part of a 
triple system with ~76 yr orbit and a mass function of ~ O.OO1O3M . 
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2. Observations 



The observations were acquired with a High-Speed Three Channel Pho- 
tometer attached to the 48 cm Cassegrain telescope at Ege University Ob- 
servatory. The observations were continued in UBVR bands. Some basic 
parameters of the program stars are listed in Table 1. The names of the 
stars are listed in the first and second columns, while J2000 coordinates are 
listed in the third and fourth columns. The V magnitudes are in the next 
column, and the spectral types of the stars are listed in the last column. 
Although the variable and comparison stars are very close on the sky, dif- 
ferential atmospheric extinction corrections were applied. The atmospheric 
extinction coefficients were obtained from observations of the comparison 
stars on each night. Moreover, the comparison stars were observed with the 
standard stars in their vicinity and the reduced differential magnitudes, in the 
sense variable minus comparis on, were trans formed to the standard system 



using procedures out-li ned by 



from the catalogues of 



iardie 



Landolt 



1983 



1962). The standard stars are used 



1992j). Furthermore, the dereddened 



colours of the system were computed. Heliocentric corrections were also ap- 
plied to the times of the observations. To compute the standard deviations 
of observations, we used the standard deviations of the reduced differential 
magnitudes in the sense comparisons minus check stars for each night. In 
Figure 1, the light and colour curves of KRCyg are shown. 

3. Orbital Period Analysis 

In order to study of the orbital period of KR Cyg, all the available times 
of light minimum were collected from the literature. We also observed the 



Table 1: The data for variable, comparison, and check stars. 



Star 


ID 


Alpha (J2000) 


Delta (J2000) 


my 


Sp.T. 






m 


(° ' ") 


(mag) 




KRCyg 


HD 333645 


20 09 05.6 


+30 33 01.3 


9.19 


A1V 


Comp. Star 


HD 191398 


20 08 39.4 


+30 20 15.2 


9.01 


AOV 


Check Star 


HD 333664 


20 09 18.1 


+30 13 39.4 


9.58 


AO 
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Figure 1: The V band light curve (upper panel) with both B-V (middle panel) and V-R 
(bottom panel) colour curves of KR Cyg. 
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system for two nights in 2010 and one night in 2011. Using these observa- 
tions, three times of primary minimum were obtained. Finally, 156 visual, 19 
photographic and 129 photoelectric times of mid-eclipse of KR Cyg were used 
to period study. New timings are listed in Table 2. The ot her minima times 



can b e taken from the modern database of O — C Gateway (IPaschke fc Brat 



20061 ) . The O — C residuals of KR Cyg were calculated using the following 



linear ephemeris; 



JD (Hel.) = 2455036.5239(8) + 0^.8451526(1) x E. (1) 

Then, the (O — C)j (upper panel) and (O-C)jj (bottom panel) residuals 
were plotted against the epoch number (E) in Figure 2 where the times of 
photoelectric minima are marked by filled circle. The (O — C)n residuals 
shows that the orbital period of KRCyg is changing in a long-term sinu- 
soidal form, although the visual and photographic data display large scatter 
throughout the diagram. Assuming that the most likely cause of the cyclic 
variation in the (O — C)n residuals could be light-travel-time effect (LTTE) 
due to an unseen third body in the system. The following LTTE equation, 
as given by Irwin ( jl959 ). was fitted to the (O — C)u residuals of the eclipse 



timings of the system: 



a 12 sini' 1 - e' 2 . , , , , 
At = \ sin( v + lo + e coscu \ 2 

where At is the time advance/delay due to LTTE, c is the speed of light, 
and ai2, i', e' and oj' are the semi-major axis, inclination, eccentricity and 
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longitude of the periastron of the absolute orbit of the centre of mass of the 
eclipsing binary around the center-of-mass of the triple system, respectively. 
v 1 is the true anomaly of the position of the eclipsing binary's mass-centre 
on this orbit. 

A weighted least-squares solution for two parameters (T and P) of the 
linear ephemeris of KRCyg and five parameters (ciusini', e', u', T' and P') of 
the LTTE are presented in Table 3. The observational points and theoretical 
best fit are plotted against the epoch number in the bottom panel of Figure 2. 
Assuming that the orbit of the presumed third body is circular, we obtained 
the mass function as f(M 3 ) = 0.00103M Q for the third body, using the 
Equation: 



,/ n 47r 2 . . ...o (M 3 sini') 3 

where M 1; M 2 and M 3 are the masses of the binary's components, and the 
third body, respectively. The mass of the third component can be estimated 
from Equation (3), depending on the orbital inclination. For example, the 
mass Ms^min was estimated to be 0.33 M Q for i' = 90°. Here, the total mass 
of the eclipsing system was taken as 5.51 M & . We calculated the third body 
masses at different inclination, i' values, which are shown in Table 3. The 
parameters given in Table 3 suggest that KRCyg has an eccentric orbit 
around the centre-of-mass of the third-body system with a period of ~76 
years. 
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Figure 2: The O — C variation of KRCyg (The photoelectric minima times were marked 
by the filled circles). 
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Table 2: The new times of minimum light for KRCyg. 



No. 


HJD 


Filter 


Error Min. 


1 


2455386.4355 


B 


0.0001 I 




2455386.4356 


V 


0.0001 I 




2455386.4356 


R 


0.0001 I 


2 


2455419.3832 


B 


0.0001 I 




2455419.3831 


V 


0.0001 I 


3 


2455691.5219 


B 


0.0001 I 




2455691.5220 


V 


0.0001 I 




2455691.5220 


R 


0.0001 I 



4. Light Curve Analysis 



Photometric analysis of KRCyg w as carried out using the PHOEBE 
V.0.31a software (IPrsa fc Z witter! . l2005[ ). The method used in the PHOEBE 
V.0.31a software depen ds on the method used in the version 2003 of the 



1971 



Wilsonl . 



19901). The UBVR 



Wilson-Devinney Code (IWilson &: Devinneyl . 
light curves were analysed simultaneously assuming the "detached" and "semi- 
detached binary, secondary star fills Roche lobe" configurations. In the pro- 
cess of the computation, we initially adopted the following fixed parame- 
ters: the mean temperature of the primary component (Ti), the linear limb- 
darkening coefficien ts of Xi and x 2 for various bands, the gravity- darkening 



exp onents of 



rtLucy 



1967I ) and the bolometric albedo coefficients of A 



A2 (IRucinskil . Il969l ) . The commonly adjustable parameters employed are the 
orbital inclination (i), the mean temperature of the secondary component 
(T 2 ), the potentials of the components and f2 2 ) and the monochromatic 
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Table 3: The parameters derived from O — C analysis of KRCyg 



Parameter 



Value 



T (HJD) 

P (day) 
K (day) 
P' (year) 
e' 

J (deg) 
ai 2 sini' (AU) 
f(m) (M ) 
M 3 (z'=15 ) 
M 3 (i'=30°) 
M 3 (i'=45°) 
M 3 (i' =60°) 
M 3 (z'=75°) 
M 3 (i'=90°) 



2455036.5239±0.0008 
0.8451526±0.0000001 
0.010±0.001 
75.6±0.6 
0.51±0.03 
28.65±0.36 
1.8044±0.0041 
0.00103±0.00009 
1.42 M & 
0.68 M Q 
0.48 M 
0.38 M 
0.34 M 
0.33 M 
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luminosity of the primary component (Li). The third light (L 3 ) was used 
also as free parameter to check for the third light contribution suggested by 
the period analysis. 

We determined the colours of the system at phase 0.75 as follow: U — 
B = m 101, B — V — m .156, V - R = m 128. We derived dereddened 
colours as (U - B) Q = -0 m .400, (B - V) = -0 m .12. Then, we took JHK 
magnitudes of the s ystem (J = 9 m .27 8, H = 9 m .169, K = 9 m .098) from the 
2MASS Catalogue (ICutri et al.l . 120031 ) . Using these magnitudes, we derived 
dereddened colours as a (J - H) = m .109 and (H - K) = m 071 for 
the system. Using the calibrations given by Tokunaga (2000), we derived the 
temperature of the primary component as 11400 K and 11800 K depending on 
the UBV and JHK dereddened colours, respectively. Both of them indicate 
the same spectral type which is B8V. We adopted the mean temperature 
of the primary component as 11400 K for the light curve analyse. To find 
a reliable photometric mass ratio, the solutions are obtained for a series of 
fixed values of the mass ratio from q = 0.30 to 0.6 in increments of 0.05. The 
sum of the squared residuals, (Sres 2 ), for the corresponding mass ratios are 
plotted in Figure 3, where the lowest value of (Sres 2 ) was found at about 
q = 0.45. The photometric elements for the mass-ratio of 0.45 are listed in 
Table 4 and corresponding light curves are plotted in Figure 4 as continuous 
line. Moreover, the 3D model at phase 0.2 and Roche geometry of the system 
is displayed in Figure 5. Using the contribution of the third body to the total 
light of the system in U, B, V and R band, we calculated the instrumental 
U — B, B — V and V — R color indices of the third component as 0.027, 
—0.022 and —0.163, respectively. These color indices that the third body 
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could be accepted as a normal star. Using the Mass-Luminosity relation of 
the main sequence stars, LaM 3 ' 5 , we calculated the values of the third body 
luminosities. We compared these luminosities with the contribution of the 
third component in the light curve analysis and we obtained the mass of the 
third body as 0.48 M , corresponding to the inclination value i' = 45°. 

Although there is no available radial velocity curve of the system, we 
try to estimate the absolute parameters of the components. Cons idering its 



spectr al type, we estimate the mass of the primary component from lTokunaga 



(120001 ) as approximately 3.8 M Q , and the mass of the secondary component 
is calculated from the estimated mass ratio of the system. Using Kepler's 
third law, we calculate the semi- major axis (a), and then the mean absolute 
radii of the components. All the estimated absolute parameters are listed 
in Table 5. In Figure 6, we plot the components in the mass-ra dius plane . 



The Mass- Radius relationships for the ZAMS and TAMS model (IPols et al. 



19981 ) also shown by the line and dotted line, respectively. In the figure, 



the open circle represents the secondary component, while the filled circle 
represents the primary component. Both of the components are located in 
the main-sequence band. 



5. Summary 

We obtained new UBVR observations for KR Cyg and analysed the UBVR 
light curves to reveal a orbital parameters solution including its Roche con- 
figuration. A detailed analysis of the orbital period changes of KR Cyg was 
made using its O — C diagram constructed from all available times of eclipse 
minima. The findings from this study are summarized as follows: 
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Figure 3: The variation of the sum of weighted squared residuals versus mass ratio. 
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Table 4: The parameters obtained from the light curve analysis. 





Detached 




Semi-detached 


Parameter 




Value 




Q 




0.45 






87.8±0.1 




87.74±0.04 


Ti (K) 




11400 




T2 (K) 


6171±19 




6055±21 


Qi 


2.856±0.003 




2.872±0.002 


^2 


2.718±0.003 




2.778±0.002 


U/L T (U) 


0.961±0.004 




0.969±0.003 


U/L T (B) 


0.951±0.004 




0.961±0.003 


U/L T (V) 


0.921±0.004 




0.933±0.003 




0.891±0.004 




0.905±0.003 


l 3 /l t (c/) 


0.008±0.003 




0.005±0.002 


L 3 /L T (B) 


0.009±0.002 




0.005±0.002 


L3/MV) 


0.008±0.002 




0.005±0.002 


L 3 /L T ( J R) 


0.007±0.002 




0.004±0.002 


#2 




1.0, 0.5 








1.0, 0.5 




Xl,bol, %2,bol 




0.719, 0.643 




Xl,U, %2,U 




0.587, 0.869 




Xl,B, %2,B 




0.654, 0.828 




x l,V, X 2,V 




0.560, 0.745 




%1,R) %2,R 




0.460, 0.653 




< ri > 


0.4501±0.0006 




0.4456±0.0005 


< r 2 > 


0.3265±0.0008 


43 


0.4187 
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Figure 4: KR Cyg's light curves observed in UBVR bands and the synthetic curves derived 
from the best light curve solutions in each band. 
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Figure 5: The 3D model at phase 0.2 (upper) and Roche geometry (bottom) of KRCyg. 

Table 5: The estimated absolute paramaters derived for KRCyg. 

Parameter Primary Secondary 

Mass (M ) 3.80 1.71 

Radius (R e ) 2.99 2.17 

Luminosity (L ) 135.4 6.10 

M bot (mag) -0.59 2.78 

log g 4.10 4.00 
d (pc) 452.16 
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Figure 6: The places of the components of KRCyg in the Mass-Radius plane. The con- 
tinuous and dott ed lines represent the ZAMS and TAMS theoretical model developed by 
Pols et all |l998l ). 
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(i) New physical and geometrical parameters of components of the system 
have been derived. Orbital parameters indicate that KR Cyg is a near-contact 
binary system with a large temperature difference of approximately 5230 K. 

(ii) KR Cyg is a semi-detached binary system, in which the primary (mas- 
sive and hotter) component is slightly smaller than its Roche lobe, while the 
cool component is slightly greater than it. It exhibits a j3 Lyrae typ e light 



curve , and fulfils most of the properties of the near-contact systems (jShawl . 



1990). 

(iii) The system can be classified as near-contact binaries of FO Vir type. 
However, it must be noted that the existing observational material does not 
allow a precise classification for the system. 

(iv) O — C analysis indicates that KR Cyg is probably an interesting triple 
eclipsing system showing the slow LTTE caused by a third body orbiting with 
a period of ~76 years. Furthermore, the significant contribution of the third- 
body found in the light curve analysis. Thus, it confirms that third body 
may produce the sinusoidal variation in the O — C residuals. 

(v) In future work, spectroscopic observations should be made to obtain 
radial velocity curves, which will allow a better discussion of the absolute 
dimensions and the evolutionary status of KRCyg. New timings of this 
eclipsing binary are also necessary to improve the LTTE parameters derived 
in this paper. 
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